Heavy Meson Description with a Screened Potential 



P. Gonzalez 1 ^, A. Valcarce 1 ^, H.Garcilazo^, and J. Vijande^ 

(1) Dpto. de Fisica Teorica and IFIC 
Universidad de Valencia - CSIC, E- 46 100 Burjassot, Valencia, Spain 
(2) Grupo de Fisica Nuclear 
Universidad de Salamanca, E-37008 Salamanca, Spain 
(3) Escuela Superior de Fisica y Matemdticas 
Instituto Politecnico Nacional, Edificio 9, 
07738 Mexico D.F., Mexico 



Abstract 



We perform a quark model calculation of the bb and cc spectra from a screened 
funnel potential form suggested by unquenched lattice calculations. A con- 
nection between the lattice screening parameter and an effective gluon mass 
directly derived from QCD is established. Spin-spin energy splittings, lep- 
tonic widths and radiative decays are also examined providing a test for the 
description of the states. 
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I. INTRODUCTION 



Heavy quarkonia, bb and to a lesser extent cc (for tt the weak decay rates are too large for 
these resonances to be observed), are ideal systems to test quark potential models derived 
from QCD. Actually for the low-lying states (the ground state for cc and the ground and 
first excited states for bb) one can test the theory itself at the perturbative regime [1]. This 
may provide a quantitative estimate of the approximations involved when using a potential 
model interaction derived from QCD to describe heavy quarkonia. 

At very short distances a nonrelativistic coulombic potential, with strength proportional 
to a s , the quark-gluon-quark coupling constant, is derived perturbatively from the one- 
gluon exchange interaction in QCD. Relativistic corrections involve spin-independent, spin- 
spin, spin-orbit and tensor terms. Radiative corrections have been estimated as well. At 
long distances nonperturbative effects involving multigluon interactions, in particular those 
related to confinement, have to be considered. Lattice calculations in the so-called quenched 
approximation (only valence quarks) derive a heavy quark long range potential linearly 
dependent on the interquark distance [2]. 

The resulting interaction can then be parametrized by a funnel shape coulomb+linear 
potential. Such a potential has been applied to get a good description of the spectrum of 
spin-triplet states for cc and for the first excited states of bb [3]. When completed with some 
relativistic plus radiative corrections it can give account, through a phenomeno logical fit of 
the mass parameters, for the whole meson spectrum ranging from heavy to light mesons, as 
well as for the electromagnetic, strong and weak decays [4,5]. 

Further QCD corrections to the coulomb+linear form can be incorporated in the potential 
through unquenched (valence+sea quarks) lattice calculations [2] . The spontaneous creation 
of light quark pairs may give rise to a breakup of the color flux tube between the two heavy 
quarks. It has been proposed that this translates into a screened funnel potential [2,6,8], 
say the potential does not rise continuously with the interquark distance but it saturates at 
the splitting energy of the heavy quark pair. Although string breaking can not be said to 
be definitely confirmed through lattice calculations [7] a quite rapid cross-over from a linear 
rising to a flat potential is well established in SU(2) Yang-Mills theories [8]. The analysis of 
the consequences of such a screening in heavy quarkonia is our main goal in this article. To 
this purpose we discuss in Sec. II the screened form of the potential and connect it to an 
effective running gluon mass in QCD. In this form we can obtain from the underlying theory 
the value of the screening parameter for a given energy scale. Then, via a sensible choice 
of this scale, we calculate in Sees. Ill to V the spectra, spin-spin splittings, leptonic widths 
and radiative decays of bb and cc mesons. The dependence of our results on the chosen scale 
is studied in Sec. VI. Finally in Sec. VII we summarize our results. 

II. THE SCREENED FUNNEL POTENTIAL 

The long range behavior of the heavy quark potential can be assigned to the saturation 
property of the running QCD coupling constant, a s (Q 2 ), at the low momentum scale. In 
a particular framework it has been shown [9], through a solution of the Schwinger- Dyson 
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equation, that when decreasing Q 2 , a s does not increase continuously but instead it saturates 
at a constant value a s (0) ~ 0.8. Explicitly: 

Air 

«s(Q 2 ) = F — T (1) 

f3 \n[(Q 2 + AM 2 (Q 2 ))/A 2 ] 

where (3$ = (33 — 2n/)/3, being nf the number of flavors with mass much smaller than Q. 
A is the QCD scale parameter for which we shall take 300 MeV hereforth and M g is an 
effective gluon mass given by: 



In 



(Q 2 + 4m 2 g )/A 2 
In (4m2/A 2 ) 



-12/11 



m 2 aq 2 ) = m 2 I L ;~ J i (2) 



with m g ~ 300 MeV. Hence the effective gluon mass, M g , runs from at Q 2 — > oo to m g 
at Q 2 — > 0. Correspondingly a s runs from at Q 2 — > oo (asymptotic freedom) to ~ 0.8 at 
Q 2 — > 0. By implementing the one-gluon-exchange diagram with a s (Q 2 ) as given above but 
substituting M 2 (Q 2 ) by M|(Q 2 -> 0) ~ A 2 , a long range (r >> A *) linear potential comes 
out [1] (one should keep in mind that fine structure splittings might require a scalar [10] or 
a mixed scalar- vector structure for the confinement [11]). 

This linear tendency gets modified by qq pair creation in between the heavy quarks 
since a screening of the quark color charge at large distances takes place. Moreover the 
effective coulomb coupling is also affected by the presence of sea quarks. These effects were 
parametrized in exploratory lattice studies by a screened funnel potential [6]. For the sake 
of simplicity we shall follow this parametrization though it does not reproduce the rapid 
turnover around 1 fm from linearly rising to flat potential suggested by modern lattice 
results [8]. For a meson the potential has the form: 

™-(-S)(nF) 

where \x is the screening parameter (with units of an inverse length) and the hat over the 
parameters in the potential distinguishes them from the nonscreened case: V(r) = or - 
V(r) behaves like a coulomb potential for r — > whereas it tends to a//i for r — > oo. 

Regarding the value of \x let us realize that for /i — > one formally recovers the non- 
screened potential with the long distance linear behavior one had when doing M g (Q 2 ) ~ A. 
On the other hand the confining part of the potential can be written in the form 

— a e~ Mr _ _ e~^ r 

V conf (r) = or = V confer -> oo) - or (4) 

/i fir fir 

what suggests that fi should be connected to an exchanged (Yukawa type) mass. From these 
considerations we propose this mass to be the effective gluon mass through the identification 

fi = A — M g (5) 

so that for a given scale specified by Q 2 



4 o, 
3 r 



3 



-jxr 



V C onf(r) = -- or- — 
H fir 



= or 



e -(A-M 9 (Q 2 ))r 



(6) 



This identification establishes a deep connection between the saturation of the coupling 
constant and the interquark pair creation mechanism both effects governed by M g {Q 2 ). 
Therefore // runs with Q 2 so that ~ /i(Q 2 = 0) < /i(Q 2 ) < /i(Q 2 -> oo) ~ 1.52 fin -1 . 

The splitting energy of quark and antiquark in the meson is given by o/ /i, the maximum 
value of the potential. No bound state can be found for higher energies due to the breaking 
of the color flux tube and the most favored subsequent decay in multimeson states. We 
shall assume that this splitting energy does not depend on flavor or, more precisely, that 
it does not depend on the scale Q 2 . Hence o should also run with the scale Q 2 in such 
a way that o(Q 2 )/fi(Q 2 ) ~ constant. For us o will be a free parameter to be determined 
phenomenologically at a chosen scale. As a value for comparison we could think of the 
string tension value for the nonscreened case, o ~ 1000 MeV fur 1 [12]. A more physical 
comparison can be done by evaluating, from o, /i, and a s , a hadronic scale R defined 
through the force between static quarks at intermediate distances as F(R ) Rq = 1.65. For 
the nonscreened case Rq ~ 0.5 fm [13]. Concerning the other parameter of the model, a s , 
since we have established a link between saturation and screening we shall assume a s ~ a s 
as given by Eq. (1) at the chosen scale. 

Our meson hamiltonian will then be written, up to an additive constant, as: 
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2m n 



2m-q 
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where the constituent quark and antiquark masses, m q and m-g, are taken as free parameters. 



III. THE bb, cc AND cb SPECTRA. 

In order to apply the screened potential to obtain the bb and cc spectra we should first 
identify the Q 2 scale in each case. For the sake of simplicity and universality let us investigate 
the possibility to describe bb and cc mesons with the same set of potential parameter values, 
say by using the same Q 2 scale. A sensible guess may be to choose Q 2 as some intermediate 
value between the square quark masses m 2 and m 2 . But it turns out that our quark model 
masses are parameters to be determined from the spectra. Nonetheless we can consider as 
a reference the central quark mass values (MS scheme) quoted by the Particle Data Book 
[14], m b = 4.2 GeV and m c = 1.2 GeV, and select Q 2 = (2.7 GeV) 2 . 

Once Q 2 has been chosen a s and \x are determined by Eqs. (1),(2) and (5). Then we 
should establish the protocol to fix the other parameters. We shall fix the quark masses and 
o to get the experimental mass for the ground states of bb and cc and the first excited state 
of bb. Indeed this corresponds to a possible definition of the mass for the confined quarks. 
Let us note that we do not expect to get the reference MS quark masses since they do not 
correspond to the same definition but quite different values will not be compatible with the 
self-consistency of our scheme. The values of the parameters used are compiled in Table I. 
The resulting spectra, square mean root radii and square quark velocities appear in Tables 
II, III and IV. 
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Let us note that our hamiltonian is spin independent what means that S = and 
5=1 states are degenerated. Assuming that the energy difference between spin singlet and 
triplet states comes out mainly due to the spin-spin interaction and taking into account the 
magnitude of these splittings and the fact that the matrix element for the spin-spin operator 

— * — * 

Si ■ Sj is three times bigger for singlets than for triplets, we shall consider our quark model 
states as close description to spin triplets. 

The good agreement of the mass prediction with data for bb and cc up to 1 GeV excitation 
energy allows an unambiguous identification of experimental J p = 1~ states: Is, 2s, Id, 3s 
and 4s for bb, and Is and 2s for cc. This is reinforced by the calculation of leptonic width 
ratios, [T e+e - /T e+e -(ls)], (see Sec. IV). The same analysis suggests that for cc our Id and 
3s states must have some mixing in order to explain the experimental values of the leptonic 
width ratios, bigger for Id and smaller for 3s than our approximated results. 

For p states and higher energies the situation is much less clear and the interpretation 
of our results requires some additional considerations. First, relativistic effects, specially for 
the cc system (see the (v 2 /c 2 ) column in Tables II, III and IV) may be quite relevant: higher 
order kinetic energy terms increase in importance when increasing the excitation energy, the 
spin-orbit force can be responsible for the p splittings and for differences between L = and 
L = 2 states, the tensor interaction can induce s — d mixing, the color magnetic potential 
breaks the degeneracy with 0~ states. Second, coupling to D and B mesons may become 
relevant [15]. Third, there can be some bias in the results associated to the Q 2 scale assumed. 
Anyhow a glance at the experimental leptonic width ratios and their comparison to our 
calculation, Tables V, suggests that T(10860) and T(11020) should be mostly s states since 
for d states the predicted ratios would be three orders of magnitude smaller. For ■0(4160) 
and -0(4415) an s state association seems also to be favored, Tables VI and IX. 

Concerning cb the commented results and the good prediction for the ground state gives 
us also confidence about the predicted masses of its lower excitations. 

The quality of the spectra we get is quite similar to that of other quark model calculations 
that use strict (nonscreened) confinement [3,4]. The major difference concerns the possible 
number of quark- ant iquark bound states which is finite in our case and infinite for the strict 
confinement models. This also affects the pattern of energy differences but mostly for higher 
excited states. Unfortunately the uncertainties associated to the calculation of these states 
are such that a comparison to data cannot discriminate against any model. The upper 
energy cut for the spectrum is given in our model by m q + trig + a/ fi. From the values of 
the parameters we obtain a splitting energy a/fj, ~ 2070 MeV and the energy thresholds 
E th {bb) ~ 11400 MeV, E th (cb) ~ 7973 MeV and E th (cc) ~ 4547 MeV. Although these values 
have to be taken very cautiously (actually they are significantly lower than the ones obtained 
from a Regge analysis [17]) since they depend on the somewhat arbitrary Q 2 scale (see Sec. 
VI) the qualitative consequence is clear: the mere existence of a threshold means that a 
detected state above it would correspond to a wide resonance, a glueball or an exotic state. 
Moreover, it could even correspond to a narrow resonance provided the existence of selection 
rules preventing its strong decay. 

It is worth to remark that our quark mass values are smaller than the ones obtained 
from strict confinement models (mj ~ 5000 MeV, m c ~ 1700 MeV) may be indicating that 
relativistic effects associated to quark pair creation can be taken into account in an effective 
way, following the philosophy of quark model calculations, via a bigger nonrelativistic quark 
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mass. Concerning the value of a (1470 MeV fm" 1 ) it is quite big as compared to commonly 
accepted values for the string tension a. However the hadronic scale provided by our model 
Rq ~ 0.43 fm is much closer to the nonscreened case. We should also keep in mind that a is 
an effective parameter which may give account of some nonconsidered corrections and that 
a runs with Q 2 in such a way that the lower the Q 2 the lower the a (see Sec. VI). 



A. Singlet-triplet states 

As the screened potential we use has no spin-dependence, spin singlet (S = 0) and triplet 
(S — 1) states are degenerated. If we consider the color- magnetic spin-spin term coming 
from the one-gluon-exchange diagram in QCD and screen it in the same manner than the 
rest of the potential, i.e., 

4 87r - - f 1 - e~^ r \ 

V ss = S t ■ S 3 6(f) e — (8) 

3 3m q mg \ fir J 

we realize that when using it in the Schrodinger equation the presence of 8(f) gives rise 
to an unbound low-energy spectrum. This difficulty can be overcome substituting the 8(f) 
by a spreading function depending on a characteristic length. However as this substitution 
represents also some kind of screening, avoiding the contact interaction, the inverse screening 
length, /i, appearing in Eq. (8) could not coincide with the value of \i for the rest of the 
potential. Then we would be left with two new parameters, the spreading length and a 
modified inverse screening length, to fit one established and two to be confirmed experimental 
splitting energies: M(J/i/j) - M[rj c (ls)], M[^(2s)\ - M[r] c (2s)] and M[T(ls)] - M[rj b (ls)], 
what does not seem to be a big deal. Instead we can try to preserve the simplicity of our 
model and apply first order perturbation theory. Thus the splitting energies obtained from 
J/ip, ^>(2s) and T(ls) are expressed as 

where \1/ (0) stands for the wave function at the origin for J/ip, ip(2s) and T(ls) respectively. 
So we predict: 

M(J/iP) - M[?7 c (ls)] ~ 118.6 MeV (10) 



M[T(ls)] - M[r] b (ls)] ~ 104.3 MeV (11) 

M[ip(2s)} - M[r] c (2s)] ~ 69.4 MeV (12) 
to be compared to the experimental data [14]: 

M(J/iP) - M[?7 c (ls)] = 117.9 ± 2.1 MeV (13) 

M[T(ls)] - M[r] b (ls)} = 160 ± 40 MeV (14) 
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and [18] 



M[ip(2s)} - M[r] c (2s)} = 64 ± 12 MeV . 



(15) 



Though we get a remarkable agreement, as the magnitudes of the splittings are not much 
smaller than the nonperturbed Schrodinger equation eigenvalues we have [620 MeV for J/i/j, 
1202 MeV for ip{2s) and 131 MeV for Y(ls)], specially for bb, the first order perturbed 
values should not be taken for granted. Anyhow we can tentatively make some additional 
predictions for the next two splittings: 



M[T(2s)] - M[q h {2s)} ~ 50 MeV 
M[-0(3s)] - M[n c {3s)} ~ 52 MeV 



(16) 



(17) 



IV. LEPTONIC WIDTHS 



A complete calculation of V (vector meson) — > e + e~ widths involves radiative and rel- 
ativistic contributions out of the scope of our quark model estimation. Fortunately these 
corrections can be factorized [19]. Thus for s states we can write the leptonic width r e + e - 
as 



r e+e -(ns) = rfl_(ns) 



16a s 
Ik 7 



+ A(ns) 



(18) 



The second ( — 16o; s /37r) and third [(A(ns)] terms in the parenthesis on the right hand 
side stand for the leading order radiative and higher order radiative+relativistic corrections 
respectively and 



r 



(0) 



M = — r^|^(0)| 2 



M 2 

ns 



(19) 



where e q is the electric quark charge, a the fine structure constant, M ns the mass of the ns 
state and \l/ riS (0) its wave function at the origin. 

From Eqs. (9), (18), and (19) we can establish an approximate relation between the 
correction factors and the quark masses through experimental quantities: 



m 2 c 


i 16a s 
. 3tt 


+ A(ls)~ 


cc 


[r e +e- 




[M(T(la)) 


- M( Vb (ls))} 


r m(j/^) i 2 


m 2 h 


1 16a s 
3tt 


+ A(ls)" 


bb 


[r e + e - 




[M( J/^) - 


-M( Vc (ls))] 


M(T(ls))_ 



(20) 



From Eqs. (18) and (19) it is also clear that a comparison of the calculated 1^+ (ns) 
to data will give us the effective magnitude of A needed in our model. So for bb we obtain 
rf + } e _(ls) = 2.13 keV and A(ls) = 0.16 and for cc we have r^ e _(ls) = 6.37 keV and 
A(ls) = 0.36. We observe that higher order corrections are much more important for the 
cc system as expected. 
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As A includes relativistic corrections it depends on the ns state. Inasmuch this depen- 
dence is not very strong we could assume a mean value for all these states. Then we could 
consider the ratio: 



Y e+e -(ns) 



[ns 



\Vns(0)\ 2 Ml 



r e+e -(is) r ° + Ui s ) |*i a (o)| 2 mi 



(21) 



where correction factors cancel. Indeed as far as the masses of the states agree with data 
this is a test of the ratio of the wave functions at the origin. 
Moreover from Eqs. (9) and (21) we can predict 





-e-(2s)" 






-e-(la). 


CC 



M(r) c (2s)) 



M(J/^) - M(n c (ls)) 



close to the experimental data, 



(2s) 



Is) 



0.41 ±0.07, 



M{j){2s)) - M(n c (2s)) 
M(J/4>) - M( Vc (ls)) 

For d states we use: 

r!tH = 



M(J/^) 



1 2 



M{i){2s)) 



0.55 ±0.04. 



2he 2 q a 2 



KM 



(22) 

(23) 
(24) 

(25) 



where -R^(O) stands for the second derivative of the radial wave function at the origin. 

Our results, shown in Tables V and VI, are quite encouraging, they seem to confirm the 
correctness of our assumptions about the validity of the prediction of ratios and support the 
state identification done in Sec. Ill (let us comment that there could be some bias due to 
the Q 2 scale value chosen, see Tables VIII and IX, Sec. VI). 



V. El AND Ml DECAY WIDTHS 

According to the identification of our states with spin triplets done in Sec. Ill, we 
shall center in El decays in the bb and cc sectors involving only triplet states, for instance 
Xcj(lp) "~ iJ/^i T(2s) — > 7X6j(lp), .... Then for i — > 7/, [i(f) stands for the initial(final) 
meson] , using a single quark current operator and the nonrelativistic approximation, we can 
write for the width 

If/ = ^e 2 q akf f (2J f + l)D 2 f (26) 

where fcj/ is the photon energy or momentum, Jf is the total angular momentum of the final 
meson and Dif the transition matrix element 

f°° 3 
D if = / drui(r) — 

JO Kif 



-Ji 



k if r 



u f (r) 



(27) 



8 



being Uij(r) the reduced radial wave functions for the initial and final mesons respectively 
and jo an d ji spherical Bessel functions. In the limit k^r — > one recovers the dipole form: 

D if -> / drui(r)ru f (r) (28) 
J o 

In our simple quark model all the J states for p waves, \Ji are degenerated. Let us 
assume that this degeneration can be dropped out by the consideration of some perturbative 
interaction such as a spin-orbit one and that first order perturbation theory can give us the 
right masses. This is quite reasonable since the needed splittings are about one tenth of the 
energy eigenvalues. Hence we shall use our wave functions (first order perturbation theory 
does not affect the wave functions) altogether with the experimental masses to get the El 
rates. 

Certainly there are relativistic and radiative corrections to the formulas above (involving 
wave function corrections also). We shall take as a criterium to estimate its importance the 
value of kifR (R can be assimilated to the mean radius of the initial meson) as compared to 1. 
Then for the transitions we consider next and from Tables II and III we have (kifR) b r ~ 0.35 
and (kifR) - > 1. Hence we expect a much more reasonable description for bb than for cc. 
This expectation is confirmed by our results as can be checked in Table VII: the bb data are 
quite well described whereas the cc widths deviate by a factor of about 2. 

The much bigger quality of the results we obtain for bb as compared to the ones obtained 
with a simple nonscreened linear potential [3] points out the importance of screening effects. 
On the other hand by taking into account that with a more refined model without screening 
[4] both bb and cc widths are correctly described we can interpret that the effect of screening 
is taken into account in this model in an effective manner through the fitted parameters. 

Less satisfactory from the experimental point of view is the situation for Ml decays since 
there are only a few experimental cc data available [J ftp — > ^rjdls) and ip(2s) — > •yrj^ls)}. 
The nonrelativistic theoretical expression for the width is: 

r " 1 = y(2^) 2/i ^ (2J ^ + 1)M ^ (29) 

with the transition matrix element Mif given by 




and where we have used an effective magnetic factor to take into account corrections to 
the quark magnetic moment beyond the Dirac particle term. 

As the J/ip — i] c (ls) and the ip(2s) — i] c (2s) mass splittings are reasonably well reproduced 
by first order perturbation theory we follow the same philosophy as above and assume the 
same wave function for J/tp and r) c (ls) and analogously for ip(2s) and rj c (2s). Let us realize 
that if we fit y? q from data this factor can effectively incorporate corrections of the same 
type mentioned for the El case as well as corrections that would come from a better wave 
function description (specially for singlet states). 

Then, for J/ip — r) c (ls), by using that M^f ~ 1 we can estimate Y j/^ 1Vc ^ ~ 4.26 p%.. 
By comparing with the experimental value, Tj/^ 1Vc ^ s ) — 1.2 ± 0.4 keV we get fj, c ~ 0.52. 
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We can compare this value with the corresponding to the second order correction to the 
magnetic moment operator [20] (1 — 2v 2 /3c 2 ) ~ 0.8. Since our /i c has a very effective 
character the value determined from i/j(2s) — > 7^ (Is) could be different. If keeping in mind 
these caveats we assume the same /x 2 we can predict r^( 2s )^ 7??c (i s ) — 1.24 keV to be compared 
to r^( 2s )- >7 ^ c (i s ) = 0.8 ± 0.3 keV. 

Let us also realize that we can get rid of the presence of \x q by taking the ratio of the 
widths. Thus we can predict for the non-measured yet ip(2s) — > i] c (2s) transition: 

r^(2s)^7?7 e (2s) ^ / k^(2s)^r) c {2s) 
^J/il>-+yr) c (ls) \ ^J/V^7»?cO) 

since Mif ~ 1 for both transitions. 

Finally for bb, since we expect relativistic corrections to be much less important, by 
assuming pL b ~ (1 — 2t> 2 /3c 2 ) = 0.94 (very close to 1, the value for a Dirac particle), we 
predict r T (i s )^ 7 f fc (i s ) = 0.17 keV (we have used M[rj b (ls)] = 9300 MeV [14]; by including 
the experimental range M[r] b (ls)] = 9300 ±20 ±20 MeV we have r T(ls) ^ m(ls) = 0.07-0.33 
keV). 

VI. THE Q 2 SCALE 

Although our results are good enough to provide an ad hoc justification for the universal 
Q 2 chosen for bb and cc there can be convenient to try to justify it from more physical 
grounds. If we write Q 2 = —q 2 , being q is the momentum transfer between quark and anti- 
quark, we have, in the center of mass system, Q 2 = 4p 2 where p q is the quark trimomentum 

p 2 = m 2 v 2 / (l — v 2 ^j. By taking as a reference the values of quark masses and velocities 

previously obtained for Is states we get Q 2 ^ ~ (2.9 GeV) 2 and Q 2 cZ ~ (2.4 GeV) 2 , pretty 
close to the Q 2 scale value we guessed Q 2 = (2.7 GeV) 2 . Let us also note that if we had used 
the nonrelativistic expressions for the trimomentum we would have got Q 2 -^ ~ (2.8 GeV) 2 

and — (1-3 GeV) 2 what reflects the much more relativistic character of cc. Anyhow 
as we are using a nonrelativistic scheme it can be of help to analyze the variation of our 
results when a Q 2 ~ (1.3 GeV) 2 scale is chosen for cc. Actually it is possible to choose 
Q2_ = ml = (1357 MeV) 2 , so that pi = 0.556 fin" 1 and a = 1175 MeV fin -1 , very close to 
the string tension value. Furthermore af/i — 2113 MeV, very close to the value obtained 
with the former universal scale what seems to confirm our assumption about its constancy. 
The results are shown in Tables VIII and IX. Due to the bigger value of the quark mass 
the upper energy threshold increases with respect to our former value (4827 MeV vs 4547 
MeV). From the comparison with Tables II, III, IV, V and VI, it is clear that the effect of 
lowering the scale translates in a smaller size and a bigger wave function value (for ns waves) 
at the origin. The rates T e + e -(ns)/T e + e -(ls) get reduced approaching data. El decay rates 
change also, although by a small amount, in the direction pointed by data. These tendencies 
though not very significant seem to favor for cc a lower scale than the one used in Sec. III. 
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VII. SUMMARY 



We have studied heavy quarkonia from a simple screened funnel quark potential model. 
The form of the potential is suggested by unquenched lattice calculations. The screening 
inverse length, i. e. the parameter of the potential controlling the onset of screening ef- 
fects, has been related to an effective gluon mass derived from QCD. As such a mass is also 
responsible for the saturation property of the running coupling constant in QCD an under- 
lying connection between saturation and screening comes out. The model provides a quite 
accurate description of bb masses, spin-spin splittings, leptonic widths and radiative decays. 
This may be pointing out the relevance of string breaking effects although the more genuine 
predictions associated to screening, say the finiteness of the quark-antiquark bound state 
spectrum, its corresponding upper energy threshold and the breaking color flux tube energy 
between quarks, remain to be tested in the future since current high energy excitation data 
do not allow to extract neither quantitative nor qualitative conclusions about their validity. 

For cc the implementation of other relativistic corrections, apart from screening, seems 
to be essential to approach data as confirmed by the accurate description provided for more 
refined models. There could be, in our opinion, interesting to implement screening effects 
in such kind of models to try to disentangle the role played by relativistic effects of different 
character. 

Let us say finally that the model can be also applied to light hadrons, mesons and 
baryons, but at the price of introducing more effectiveness in the values of the parameters, 
losing the stringent connection established with the underlying theory. 

ACKNOWLEDGMENTS 

We are specially grateful to Dr. M.A. Sanchis for providing us with phenomenological 
information. This work has been partially funded by Ministerio de Ciencia y Tecnologia 
under Contract No. BFM2001-3563, by Junta de Castilla y Leon under Contract No. SA- 
109/01, and by EC-RTN (Network ESOP) under Contract No. HPRN-CT-2000-00130. 



11 



REFERENCES 



[1] F.J. Yndurain, hep-ph/9910399. 
[2] G.S. Bali, Phys. Rep. 343, 1 (2001). 

[3] E. Eichten, K. Gottfried, T. Kinoshita, K.D. Lane, and T.M. Yan, Phys. Rev. D 17, 
3090 (1978); ibid 21, 203 (1980). E. Eichten and C.Quigg, Phys. Rev. D 52, 1726 (1995). 
[4] S. Godfrey and N. Isgur, Phys. Rev. D 32, 189 (1985). 

[5] R.K. Bhaduri, L.E. Cohler, and Y. Nogami, Nuovo Cimento A 65, 376 (1981). 
[6] K.D. Born, E. Laermann, N. Pirch, T.F. Walsh, and P.M. Zerwas, Phys. Rev. D 40, 
1653 (1989). 

[7] G.S. Bali et al, Phys. Rev. D 62, 054503 (2000); B. Bolder et al, Phys. Rev. D 63, 

074504 (2001); P. Pennanen and C. Michael, hep-lat/0001015; C.W. Bernard et al, 

Phys. Rev. D 64, 074509 (2001). 
[8] PW. Stephenson, Nucl. Phys. B 550, 427 (1999). F. Knechtli and R. Sommer, Nucl. 

Phys. B 590, 309 (2000). O. Philipsen and H. Wittig, Phys. Lett. B 451, 146 (1999). 

P. de Forcrand and O. Philipsen, Phys. Lett. B 475, 280 (2000). 
[9] A.C. Aguilar, A. Mihara and A. A. Natale, hep-ph/0109223; J.M. Cornwall and J. Pa- 

pavassiliou, Phys. Rev. D 40, 3474 (1989); ibid 44, 1285 (1991). 
[10] H.J. Snitzer, Phys. Lett. B 76, 461 (1978). 

[11] Y.S. Kalashnikova, A.V. Nefediev, Phys. Lett. B 414, 149 (1997). 

[12] G. Bali, K. Schilling and A. Wachter, Phys. Rev. D 56, 2566 (1997); hep-ph/9611226. 

[13] R. Sommer, Nucl. Phys. B 411, 839 (1994). 

[14] K. Hagiwara et al, Phys. Rev. D 66, 010001 (2002). 

[15] E. Eichten, Phys. Rev. D 22, 1819 (1980). 

[16] S.E. Csorna et al, CLEO collaboration, hep-ex/0207060. 

[17] M.M. Brisudova, L. Burakovsky, and T. Goldman, Phys. Rev. D 61, 054013 (2000). 
[18] K. Abe et al, Phys. Rev. Lett. 89, 142001 (2002); S.K. Choi et al, Phys. Rev. Lett. 

89, 102001 (2002). 
[19] W. Buchmuiler and S.H.H. Tye, Phys. Rev. D 24, 132 (1981). 
[20] T.A. Lahde, C.J. Nyfalt, and D.O. Riska, Nucl. Phys. A 645, 645 (1999). 



12 



TABLES 



TABLE I. Parameters of the model. 



a s 0.317 

H (fm- 1 ) 0.71 

a (MeV fm" 1 ) 1470 

m b (MeV) 4664.5 

m c (MeV) 1238.5 



TABLE II. bb spin-triplet bound state masses (in MeV) and properties. We denote by a '*' the 
states used to fit the parameters and by a 'f experimental states with different possibilities for the 
orbital angular momentum, it could be L = and/or L = 2. Experimental data are taken from 
Ref. [14] except for that quoted by a '**' taken from Ref. [16]. p states are taken as the centroid 
of the npo, np\ and np2 experimental data. 





Mass 


Exp. 


(v 2 /c 2 ) 


(r 2 ) 1/2 (fm) 


Is 


9460* 


9460.30 ± 0.26 


0.09 


0.23 


2s 


10023* 


10023.26 ± 0.31 


0.08 


0.51 


Id 


10159 


10162.2 ± 1.6** 


0.07 


0.55 


3s 


10350 


10355.2 ± 0.5 


0.08 


0.78 


2d 


10436 




0.08 


0.82 


4s 


10588 


10580.0 ± 3.5 


0.09 


1.05 


3d 


10649 




0.08 


1.10 


5s 


10773 


10865 ±8 + 


0.09 


1.33 


Ad 


10819 




0.08 


1.40 


6s 


10921 


11019 ±8t 


0.08 


1.60 


5d 


10956 




0.08 


1.70 


7s 


11040 




0.08 


1.90 


lp 


9909 


9900.1 ±0.5 


0.07 


0.40 


2p 


10262 


10260.0 ± 0.5 


0.08 


0.70 
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TABLE III. cc spin-triplet bound state masses (in MeV) and properties. We denote by a '*' 
the states used to fit the parameters and by a 'f ' the experimental states with different possibilities 
for the orbital angular momentum, it could be L = and/or L = 2. Experimental data are taken 
from Ref. [14]. p states are taken as the centroid of the npo, np\ and np2 experimental data. 





Mass 


Exp. 


(v 2 /c 2 ) 


(r 2 ) 1 ' 2 (fm) 


Is 


3097* 


3096.87 ± 0.04 


0.28 


0.46 


2s 


3679 


3685.96 ± 0.09 


0.32 


0.97 


Id 


3794 


3769.9 ± 2.5 1 


0.31 


1.00 




4093 


4040 + i o t 




1 50 


2d 


4093 




0.31 


1.60 


4s 


4248 


4159 ± 20 1 


0.29 


2.10 


OU 


4294 




0.27 


2.20 




4397 


4415 ±6 + 


0.23 


2.95 


Ad 


4427 




0.22 


3.10 


6s 


4494 




0.21 


3.70 


lp 


3517 


3525.3 ±0.2 


0.12 


0.75 


2p 


3914 




0.20 


1.3 




TABLE IV. be spin-triplet bound state masses 


(in MeV) and properties. Experimental data 


3X6 


taken from Ref. [14]. 










Mass 


Exp. 


(v 2 /c 2 ) 


<r 2 ) 1/2 (fm) 


Is 


6362 


6400 ± 400 


0.18 


0.37 


2s 


6927 




0.20 


0.79 


Id 


7045 




0.19 


0.82 


3s 


7269 




0.20 


1.20 


2d 


7342 




0.20 


1.25 


4s 


7506 




0.20 


1.70 


3d 


7556 




0.19 


1.71 


5s 


7677 




0.17 


2.16 


Ad 


7713 




0.17 


2.24 


6s 


7799 




0.15 


2.78 


5d 


7825 




0.15 


2.89 


7s 


7884 




0.12 


3.56 


lp 


6779 




0.18 


0.60 


2p 


7163 




0.20 


1.02 
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TABLE V. Leptonic widths T e + e - (in keV) for bb. We denote by a 'f the experimental states 
with different possibilities for the orbital angular momentum, it could be L = and/or L = 2. 
Experimental data are taken from Ref. [14]. 

rf+e- (i-ff-) OVe-W rg _/rg> e -(i*) [r e+e -/r e+e -(i«)U 

Is 0.98 1.32 ±0.05 1 1 

2s 0.41 0.520 ±0.032 0.42 0.41 ± 0.07 

Id 3.7 x 10" 4 3.8 x 10" 4 

3s 0.27 seen 0.27 seen 

2d 5.8 x 10- 4 5.9 x 10~ 4 

4s 0.20 0.248 ±0.031 0.21 0.19 ±0.03 

3d 6.7 x 10~ 4 6.8 x 10~ 4 

5s 0.16 0.31±0.07t 0.16 0.24±0.06 t 

Ad 7.9 x 10" 4 8.1 x 10~ 4 

6s 0.12 0.130 ± 0.030 1 0.13 0.10±0.03 t 

5d 7.1 x 10" 4 7.3 x 10" 4 

7s 0.10 0.10 



TABLE VI. Leptonic widths T e + e - (in keV) for cc. We denote by a 'f the experimental states 
with different possibilities for the orbital angular momentum, it could be L = and/or L = 2. 
Experimental data are taken from Ref. [14]. 

rffe- OVe-W rff e _/rff) e .(i*) [r e+e -/r e+e -(is)] exp 



Is 2.94 5.3 ±0.4 1 1 

2s 1.22 2.19 ±0.15 0.42 0.41 ± 0.06 

Id 0.026 0.26±0.04t 8.8 x 10~ 3 (5.0 ± 1.1) x 10 -2 

3s 0.76 0.75 ± 0.15 1 0.26 0.15 ± 0.04 1 

2d 0.03 1.1 x 10" 2 

4s 0.43 0.77±0.23 t 0.15 0.15±0.05 t 

3d 0.03 9 x 10~ 3 

5s 0.27 0.47±0.10 t 0.09 0.09±0.03 t 

Ad 0.02 7.9 x 10~ 3 
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TABLE VII. El decay widths for bb and cc (in keV). Experimental data are taken from Ref. 



[14]. 





bb 








cc 






Transition 




r 

L exp 


Transition 




r 

- 1 - errp 


T(2s) 




1.46 


1.7 ±0.5 




- 7-Wls) 


63.94 


169 ± 51 


T(2s) 


->7X6i(l^) 


2.27 


3.0 ±0.8 


Xd(lP) 


-> 7^(1«) 


125.07 


295 ± 71 


T(2s) 




2.32 


3.1 ±0.8 


Xc 2 (lP) 


-> 7 ^(ls) 


161.64 


392 ± 73 


T(3s) 




1.71 


1.4 ±0.4 


^(2s) - 


-7X C0 (1^) 


62.56 


26 ±5 


T(3s) 


-7Xbi(2P) 


2.80 


3.0 ±0.6 


^(2s) - 


->7X C1 (1P) 


59.39 


26 ±4 


T(3s) 




3.04 


3.0 ±0.6 


V>(2s) - 


^7Xc 2 (lP) 


42.39 


21 ±4 



TABLE VIII. cc spin-triplet bound state masses (in MeV) and properties from a screened 
funnel hamiltonian with parameters: a s = 0.451, [i = 0.556 fm _1 , a = 1175 MeV fm _1 , and m c = 
1357 MeV. We denote by a '*' the states used to fit the parameters and by a 'f the experimental 
states with different possibilities for the orbital angular momentum, it could be L = and/or 
L = 2. Experimental data are taken from Ref. [14]. p states are taken as the centroid of the npo, 
npi and np2 experimental data. 





Mass 


Exp. 


(v 2 /c 2 ) 


(r 2 ) 1/2 (fm) 


Is 


3097* 


3096.87 ± 0.04 


0.28 


0.42 


2s 


3686* 


3685.96 ± 0.09 


0.29 


0.92 


Id 


3816 


3769.9 ± 2.5 1 


0.07 


0.97 


3s 


4039 


4040 ± 10 + 


0.30 


1.41 


2d 


4119 




0.15 


1.47 


4s 


4286 


4159 ± 20 1 


0.30 


1.92 


3d 


4341 




0.18 


2.00 


5s 


4467 


4415 ±6 t 


0.28 


2.50 


lp 


3542 


3525.3 ±0.2 


0.11 


0.72 
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TABLE IX. Leptonic widths T e + e - (in keV) for cc with the same model as in Table VIII. We 
denote by a 'f the experimental states with different possibilities for the orbital angular momentum, 
it could be L = and/or L = 2. Experimental data are taken from Ref. [14]. 

rffe- { l - lJ t) OVe-U rgZgg _(is) [r e+e -/r e+e -(is)] e ~ 



Is 2.30 5.3 ±0.4 1 1 

2s 0.87 2.19 ±0.15 0.38 0.41 ± 0.06 

Id 1.3 x 10~ 2 0.26 ± 0.04 1 5.5 x 10~ 3 (5.0 ± 1.1) x 10 -2 

3s 0.47 0.75±0.15 t 0.20 0.15 ± 0.04 1 

2d 1.72 x 10~ 2 7.4 x 10" 3 

4s 0.34 0.77±0.23 t 0.15 0.15±0.05 t 

3d 1.74 x 10~ 2 8.2 x 10" 3 

5s 0.24 0.47±0.10 t 0.10 0.09±0.03 t 
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